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Wo report the results of a search for a narrow resonance in electron-positron events in the invariant 
mass range of 150 — 950 GeV/c^ using 1.3 fb~^ of pp collision data at = 1-96 TeV collected by 
the CDF II detector at Fermilab. No significant evidence of such a resonance is observed and we 
interpret the results to exclude the standard model-like Z' with a mass below 923 GeV/c^ and the 
Randall- Sundrum graviton with a mass below 807 GeV/c^ for k/Mpi=0.1, both at the 95% confidence 
level. Combining with di-photon data excludes the Randall-Sundrum graviton for masses below 889 
GeV/c^ for k/Mpi=OA. 

PACS numbers: 13.85.Rm, 13.85.Qk, 14.70.Pw, 04.50.-fh 
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At hadron colliders, electron-positron pairs (ee) are a 
distinct experimental signature with a low background 
rate. Since many models introducing new physics beyond 
the standard model of particle physics (SM) predict an 
excess in ee production at a hadron collider, this channel 
has a strong discovery potential. This Letter describes 
a search for a new high mass state in ee events from pp 
collisions at ^/s = 1.96 TeV. The data used in this analy- 
sis were collected by the CDF II detector at the Fermilab 
Tevatron and correspond to an integrated luminosity of 
1.3 fb~^. The analysis also uses results from the 77 chan- 
nel described in [1[ to increase the analysis's sensitivity 
to the Randall- Sundrum (RS) graviton 

The search is optimized for new physics processes 
which produce narrow ee resonances [3|, but is other- 
wise model- independent. In addition to the above search, 
cross section times branching ratio (cr • B) upper lim- 
its Q are set for generic neutral spin-1 and spin-2 bosons. 
These a ■ B limits are then used to set lower bounds on 
the masses of specific particles predicted by new physics 
models. These particles are the Eq Z's [5| and the RS 
graviton. The Eq model unifies the forces of the SM into 
the Eq gauge group and in doing so predicts the presence 
of two additional neutral massive spin-1 bosons, referred 
to as Z's, which can mix with some arbitrary mixing 
angle. The Z'^, Z'^, Z'^, and Z'j Z's correspond to spe- 
cific values of the mixing angle and are used to bench- 
mark the model. The RS graviton is predicted by the RS 
model of warped extra dimensions which solves the hier- 
archy between the weak and Planck scales by introduc- 
ing an extra spatial dimension with negative curvature 
k. The model predicts a series of narrow neutral spin- 
2 resonances which couple to all SM particles, with the 
lowest mass resonance referred to here as the RS gravi- 
ton. The properties of this model are determined by the 
mass of the RS graviton and the ratio fc/Mp/, where Mpi 
is the reduced effective Planck scale. This ratio governs 
the couplings of the graviton to SM particles and it has 
a favored range of 0.01 to 0.1 

The most recent similar search by the D0 collaboration 
used an integrated luminosity of 260 pb"-'^ and treated the 
ee and the di-photon (77) channels as a single channel to 
perform the first dedicated RS graviton search 0]. The 
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most recent search for new physics in the ee channel by 
the CDF collaboration was a dedicated Z' search and 
used an integrated luminosity of 448 pb~^ [1]. 

This analysis is based on an integrated luminosity 
of 1.3 fb-i collected with the CDF II detector. The 
CDF II detector is a general purpose detector which is 
azimuthally and forward-backward symmetric and is de- 
scribed in detail elsewhere The relevant components 
for this analysis are the central tracking chamber (COT) 
and the central and plug calorimeters. The COT is a 96- 
layer drift chamber placed within a 1.4 Tesla magnetic 
field and is used to measure the momenta of charged 
particles within the pseudorapidity range |?7| < 1.1 [lO|. 
The COT is complemented by a silicon microstrip de- 
tector which directly surrounds the beampipe and has a 
tracking coverage of \ri\ < 2.0 [lH. The central and plug 
calorimeters are sampling calorimeters that surround the 
COT; they consist of electromagnetic (EM) and hadronic 
sections that measure the energy of particles in the range 
I77I < 1.1 and 1.2 < \r]\ < 3.6 respectively. 

The trigger used in this analysis requires two separate 
deposits of EM energy in the calorimeter and is effectively 
100% efficient for selecting ee events within the accep- 
tance of the analysis. Events are selected by requiring 
two electron candidates with Et > 25 GeV. Events are 
separated into two channels: central-central (CC) where 
both electrons are in the central EM calorimeter (CEM) 
(l^yl < 1-1) and central-plug (CP) where one electron is in 
the plug EM calorimeter (PEM) (1.2 < I77I < 3.0). Elec- 
trons in the CEM are required to have a well-measured 
track, whereas there is no tracking requirement for elec- 
tron candidates in the PEM. Due to the lack of tracking 
in the plug region there is no opposite sign requirement 
for any of the ee pairs. In the CC channel, 5% of the 
electron pairs in the signal region (Mee > 150 GeV/c^) 
are same sign, which is compatible with the fraction of 
misidentified opposite sign pairs predicted by simulation. 
Electrons are identified in an identical way to the previ- 
ously published analysis [1], with the exception that a 
photon conversion veto is applied to central electrons in 
CP events. This selection cut improves the sensitivity of 
the analysis by reducing the 77 background in this chan- 
nel. The event selection and the search method defined 
later were chosen without regard to events observed in 
the signal region to ensure a statistical robust result. 

The geometric and kinematic acceptance as a function 
of resonance mass is estimated using event samples gen- 
erated by Monte Carlo (MC) simulation. The PYTHIA 
event generator with the CTEQ5L parton distri- 
bution functions (PDF) [H] and the CDF II detector 



simulation based on GEANT 3 [lj| are used to generate 
all simulation samples unless otherwise stated. A Z' with 
the coupfings of the SM Z (SM-like Z') is used for the 
simulated spin-1 signal sample and a RS graviton with 
k/Mpi = 0.1 is used for the spin-2 sample. Both the 
Z' and RS graviton bosons are constrained to be within 
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±10% of their on-shell mass at generator level. The un- 
certainty on the acceptance resulting from the PDF pa- 
rameterization is estimated to be 2-4% using the pro- 
cedure recommended by the CTEQ Cohaboration TE^ . 
The uncertainty on the acceptance due to initial state 
radiation (ISR) is estimated to be 4% by varying the pa- 
rameters governing ISR in PYTHIA. The electron identifi- 
cation efficiency ranges from 90% to 95% per electron and 
is estimated using the simulated signal samples. These 
estimates are corrected for imperfections in the simula- 
tion by comparing the simulation with the data at the Z 
pole, resulting in a 2% systematic uncertainty. The total 
selection efficiency x acceptance over the entire signal re- 
gion is in the range of 40% to 45% for spin-1 and spin-2 
bosons. 

The most significant source of background to new 
physics in the ee channel is the SM Drell-Yan process via 
Z/j* which is an irreducible background. Jet events, such 
as di-jet or VF-|-jet events where the jets are misidenti- 
fied as electrons, represent the most significant reducible 
background. Other less significant backgrounds result 
from tt, T+T~ , WW, WZ, Wj, and 77 events. Figure [1] 
shows the invariant mass distribution for all the back- 
ground components properly normalized, together with 
the observed data for the CC and CP channels combined. 
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FIG. 1: The measured ee mass spectrum with the expected 
background for the CC and CP channels combined. The back- 
grounds are displayed cumulatively. There are no observed 
events above 550 GeV/c^. 



The SM Drell-Yan contribution is estimated using MC 
simulated events normalized to the data at the Z pole. 
By investigating the stability of the normalization factor 
a 4% systematic uncertainty is obtained on the SM Drell- 
Yan normalization. An uncertainty on the SM Drell-Yan 
shape due to PDF uncertainties is determined using the 
same method as used for the acceptance. The di-jet and 
VF-|-jet backgrounds are treated as a single background. 



referred to as the jet background. The size of this back- 
ground is estimated from a sample of jet events con- 
structed from the data identically to the signal sample 
except that at least one electron candidate is not iso- 
lated and therefore likely to be a jet. From the dis- 
tribution of the isolation of these 'electron' candidates, 
the number of jet events in the signal sample and its un- 
certainty is extracted. The shape of the jet background is 
estimated from a jet-|-F sample, where Y is either an elec- 
tron or a jet misidentified as an electron. In the region 
where the jet background is significant, the normaliza- 
tion uncertainty is the dominant uncertainty on the jet 
background. Using these methods, the jet background is 
estimated to account for 0.8 ± 0.7% and 25 ± 8% of the 
total background above 150 GeV/c^ in the CC and CP 
channels respectively. The remaining backgrounds arc all 
estimated using MC simulation normalized to the theo- 
retical next-to-leading order (NLO) or higher order cross 
section. The uncertainties in these background estimates 
are dominated by the 6% uncertainty on the luminos- 
ity ;i7]. 

A model-independent search for an excess over SM 
predictions is performed in an invariant mass range of 
150 — 950 GeV/c^. The search is optimized for a nar- 
row resonance, but still retains sensitivity to other sig- 
nals which would produce an excess over SM predictions. 
Using 1 GeV/c^ steps from M^e = 150 to 950 GeV/c^, 
the probability, referred to as the p-value, of observ- 
ing at least as many events as recorded in the real 
experiment given the expected background rate is cal- 
culated using Poisson statistics in a mass window of 
4.8 -t- 0.044 X Mee GeV/c^. This mass window is ap- 
proximately the width a narrow resonance would have 
if observed in the CDF detector, and this choice of mass 
window maximizes the sensitivity to discovering such a 
resonance for this particular analysis, as verified from 
studies of simulated events. The uncertainty on the back- 
ground estimate is treated as a nuisance parameter with 
a Gaussian distribution. The p-values for the CC and 
CP channels are combined multiplicatively and the re- 
sults are shown in Fig. [2l The minimum p-value ex- 
pected in the absence of new physics depends on the size 
of the search range, with the expected minimum p-value 
decreasing as the search range increases. The range in 
which minimum p-value is expected to occur is shown in 
Fig. 2 and is defined to include 68.3% of the minimum 
p-values centered on the median value, using one million 
simulated mass spectra. Similarly, the 3(j evidence line is 
the p-value above which the minimum p-value in 99.85% 
of the simulated mass spectra fall; any p-value lower than 
this would be taken as evidence for the presence of new 
physics. The lowest p-value observed is at 367 GeV/c^ 
and is within the expected range. It is therefore con- 
cluded that the results of this analysis are consistent with 
the SM. 

To complement the above search, a Bayesian binned 
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FIG. 2: The probability that the background alone could give 
rise to the observed number of events in a mass window equal 
to the width of a narrow resonance in the CDF detector. The 
expected range and the 3a evidence line are defined in the 
text. The region 550 — 950 GeV/c^ is not shown as no events 
are observed in data and therefore the p-value is always 1 in 
this region. 



likelihood method is used to extract limits on a ■ B{X 
ee), where the mass of X is within ±10% of its on-shell 
mass. As the acceptance of the final-state ee system is re- 
quired to extract a cross section, it is necessary to specify 
the spin of the particle. Both spin-1 and spin-2 particles 
are considered here. The likelihood is 1-dimensional with 
the signal cross section as the free parameter and the 
bin contents treated using Poisson statistics. The like- 
lihood is then convoluted with a Gaussian to allow for 
the uncertainty on the cross section from the acceptance, 
background, and luminosity estimates. The probability 
density function is formed by taking a flat prior for the 
signal cross section and is numerically integrated to ob- 
tain the limit on a ■ B{X ee). The observed limits are 
shown in Fig.[3]for the spin-1 case. The Z' model lines are 
obtained using the leading-order pythia event generator 
using the couplings in [18j, with a factor of 1.3 applied 
multiplicatively to account for NLO corrections [19| . 

For the specific case of the RS graviton, which has a 
branching ratio to 77 twice that of ee, the analysis sen- 
sitivity can be improved by combining with the 77 chan- 
nel described in jll]. The 77 events are required to have 
two photons with Et > 15 GeV, with one in the GEM 
(1 77 1 < 1.04) and the other either in the GEM or PEM 
(1.2 < \ri\ < 2.8). The GC and GP channels use inte- 
grated luminosities of 1.2 fb~^ and 1.1 fb^^ respectively, 
with the difference arising from the requirement that the 
silicon system be operational for track veto purposes for 
photons in the PEM. Selected photons have similar iso- 
lation and shower shape requirements to electrons; how- 
ever photons are required not to have an associated track. 



which ensures there is no overlap between the 77 and ee 
sample. 

The 77 channel is combined with the ee channel by 
multiplying their individual likelihoods together. The 
uncertainties on the background estimates are consid- 
ered to be uncorrelated. The uncertainties on the accep- 
tance and luminosity are taken to be 100% correlated. 
The combined limits, together with the k/Mpi vs gravi- 
ton mass exclusion region, are shown in Fig. 01 The RS 
graviton model lines are obtained using the herwig event 
generator [2^ with a factor of 1.3 applied multiplicatively 
to account for NLO corrections. 
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FIG. 3: The observed and expected limits on the a ■ B{X - 
ee) of a spin-1 particle. 



In summary a search has been made for new physics in 
the ee channel, and no significant excess over the stan- 
dard model prediction is observed. Limits are placed on 
new spin-1 and spin-2 bosons. The SM-like Z' is found to 
be excluded for masses below 923 GeV/c^ and the Z' 
bosons; the Z|, the Z'^, the Z^, and the Z'^ bosons are ex- 
cluded with masses below 729, 822, 822, and 891 GeV/c^ 
respectively. The direct limits presented here on all the 
Z' bosons surpass the corresponding indirect limits 
from LEP The RS graviton with k/Mpi = 0.1 is 

excluded for masses below 807 GeV/c^. When combined 
with the 77 channel, masses less than 889 GeV/ (? are ex- 
cluded for k/Mpi =0.1. The above limits on Z's and the 
RS graviton represent the best single-experiment direct 
limits to date. 

We thank the Fermilab staff and the technical staffs 
of the participating institutions for their vital contribu- 
tions. This work was supported by the U.S. Department 
of Energy and National Science Foundation; the Italian 
Istituto Nazionale di Fisica Nucleare; the Ministry of 
Education, Culture, Sports, Science and Technology of 
Japan; the Natural Sciences and Engineering Research 
Gouncil of Ganada; the National Science Council of the 



7 




FIG. 4: The limits on a ■ B(G ee) (top) and kjU^i (bot- 
tom) for a RS graviton in the ee, 77 channels separately and 
combined. 



Republic of China; the Swiss National Science Founda- 
tion; the A. P. Sloan Foundation; the Bundesministerium 
fill Bildung und Forschung, Germany; the Korean Sci- 
ence and Engineering Foundation and the Korean Re- 
search Foundation; the Science and Technology Facilities 
Council and the Royal Society, UK; the Institut National 
de Physique Nuclcaire ct Physique des Particules/CNRS; 
the Russian Foundation for Basic Research; the Comision 
Interministerial de Ciencia y Tecnologia, Spain; the Eu- 
ropean Community's Human Potential Programme; the 



Slovak R&D Agency; and the Academy of Finland. 



[1] T. Aaltonen et al. (CDF Collaboration), to be submitted 

to Phys. Rev. Lett. 
[2] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 

(1999). 

[3] A narrow resonance is defined as a resonance where the 
resolution of the CDF detector solely determines the ob- 
served width. 

[4] All limits presented in this Letter are at the 95% confi- 
dence level. 

[5] F. del Aguila, M. Quiros, and F. Zwirner, Nucl. Phys. 

B287, 419 (1987); J. L. Hewett and T. G. Rizzo, Phys. 

Report. 183, 193 (1989). 
[6] H. Davoudiasl, J. L. Hewett, and T. G. Rizzo, Phys. Rev. 

Lett. 84, 2080 (2000); 

B. C. AUanach et al. J. High Energy Phys. 0212, 039 
(2002). 

[7] V. M. Abazov et al. (D0 Collaboration), Phys. Rev. Lett. 

95, 091801 (2005). 
[8] A. Abulencia et al. (CDF Collaboration), Phys. Rev. 

Lett. 96, 211801 (2006). 
[9] D. Acosta et al. (CDF Collaboration), Phys. Rev. D71, 

032001 (2005). 

[10] CDF uses a cylindrical coordinate system with the z-axis 
along the proton direction, where z=0 is the center of the 
detector. Pseudorapidity defined as = — ln(tan(S/2)), 
where 6 is the polar angle while Et = i5sin(^), where 
E is the energy deposited in the cluster of calorimeter 
towers. 

[11] A. Sill et al., Nucl. lustrum. Methods A 447, 1 (2000). 

[12] T. Sjostrand et al, Comput. Phys. Commun. 135, 238 
(2001). We use pythia version 6.216. 

[13] H. L. Lai et al, Eur. Phys. J. C12 375 (2000). 

[14] S. Agostinelli et al., Nucl. lustrum. Methods Phys. Res. 
A 506, 250 (2003). 

[15] J. Pumplin et al, J. High Energy Phys. 0207, 012 (2002). 

[16] Electrons are isolated if in a surrounding cone of 0.4 ra- 
dius the Et deposited minus 1.02 x the electron Et is 
less than 3 GeV(1.6 GeV) in the GEM (PEM). 

[17] S. Klimenko, J. Konigsberg, and T.M. Liss, FERMILAB- 
FN-0741 (2003). 

[18] C. Giobanu et. al., FERMILAB-FN-0773-E (2005). 

[19] D. Acosta et al. (CDF Collaboration), Phys. Rev. Lett. 
94, 091803 (2005). 

[20] G. Corcella et al, J. High Energy Phys. 0101, 10 (2001). 

[21] W.-M. Yao et al, J. Phys. G 33, 1 (2006). 



